fileNE="result_bds.txt";

data_matrix = readmatrix(fileNE);
dx=data_matrix(:,1);
dy=data_matrix(:,2);
dz=data_matrix(:,3);
figure;
subplot(3,1,1); %定位误差
plot(dx);
title('E方向坐标');
xlabel('历元数(s)');
ylabel('E(m)');
subplot(3,1,2); 
plot(dy);
title('N方向坐标');
xlabel('历元数(s)');
ylabel('N(m)');
subplot(3,1,3); 
plot(dz);
title('U方向坐标');
xlabel('历元数(s)');
ylabel('U(m)');


dvx=data_matrix(:,4);
dvy=data_matrix(:,5);
dvz=data_matrix(:,6);
figure;
subplot(3,1,1); %速度误差
plot(dvx);
title('x方向速度');
xlabel('历元数(s)');
ylabel('m/s');
subplot(3,1,2); 
plot(dvy);
title('y方向速度');
xlabel('历元数(s)');
ylabel('m/s');
subplot(3,1,3); 
plot(dvz);
title('z方向速度');
xlabel('历元数(s)');
ylabel('m/s');

figure;
satnum=data_matrix(:,7);
pdoop=data_matrix(:,8);
subplot(2,1,1); %速度误差
plot(pdoop);
xlabel('历元数(s)');
title('pdop值');
subplot(2,1,2); %速度误差
plot(satnum);
title('卫星个数');
xlabel('历元数(s)');

%%
fileNE="result_gps.txt";

data_matrix = readmatrix(fileNE);
dx=data_matrix(:,1);
dy=data_matrix(:,2);
dz=data_matrix(:,3);
figure;
subplot(3,1,1); %定位误差
plot(dx);
title('E方向坐标');
xlabel('历元数(s)');
ylabel('E(m)');
subplot(3,1,2); 
plot(dy);
title('N方向坐标');
xlabel('历元数(s)');
ylabel('N(m)');
subplot(3,1,3); 
plot(dz);
title('U方向坐标');
xlabel('历元数(s)');
ylabel('U(m)');


dvx=data_matrix(:,4);
dvy=data_matrix(:,5);
dvz=data_matrix(:,6);
figure;
subplot(3,1,1); %速度误差
plot(dvx);
title('x方向速度');
xlabel('历元数(s)');
ylabel('m/s');
subplot(3,1,2); 
plot(dvy);
title('y方向速度');
xlabel('历元数(s)');
ylabel('m/s');
subplot(3,1,3); 
plot(dvz);
title('z方向速度');
xlabel('历元数(s)');
ylabel('m/s');

figure;
satnum=data_matrix(:,7);
pdoop=data_matrix(:,8);
subplot(2,1,1); %速度误差
plot(pdoop);
xlabel('历元数(s)');
title('pdop值');
subplot(2,1,2); %速度误差
plot(satnum);
title('卫星个数');
xlabel('历元数(s)');

%%

fileNE="result_both.txt";

data_matrix = readmatrix(fileNE);
dx=data_matrix(:,1);
dy=data_matrix(:,2);
dz=data_matrix(:,3);
figure;
subplot(3,1,1); %定位误差
plot(dx);
title('E方向坐标');
xlabel('历元数(s)');
ylabel('E(m)');
subplot(3,1,2); 
plot(dy);
title('N方向坐标');
xlabel('历元数(s)');
ylabel('N(m)');
subplot(3,1,3); 
plot(dz);
title('U方向坐标');
xlabel('历元数(s)');
ylabel('U(m)');


dvx=data_matrix(:,4);
dvy=data_matrix(:,5);
dvz=data_matrix(:,6);
figure;
subplot(3,1,1); %速度误差
plot(dvx);
title('x方向速度');
xlabel('历元数(s)');
ylabel('m/s');
subplot(3,1,2); 
plot(dvy);
title('y方向速度');
xlabel('历元数(s)');
ylabel('m/s');
subplot(3,1,3); 
plot(dvz);
title('z方向速度');
xlabel('历元数(s)');
ylabel('m/s');

figure;
satnum=data_matrix(:,7);
pdoop=data_matrix(:,8);
subplot(2,1,1); %速度误差
plot(pdoop);
xlabel('历元数(s)');
title('pdop值');
subplot(2,1,2); %速度误差
plot(satnum);
title('卫星个数');
xlabel('历元数(s)');